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, This article reviews steady-state spin densities and spin currents in materials with strong spin- 

f"^ ■ orbit interactions. These phenomena are intimately related to spin precession due to spin-orbit 

coupling which has no equivalent in the steady state of charge distributions. The focus will be 
£SJ , initially on effects originating from the band structure. In this case spin densities arise in an 

electric field because a component of each spin is conserved during precession. Spin currents arise 
because a component of each spin is continually precessing. These two phenomena are due to 
independent contributions to the steady-state density matrix, and scattering between the conserved 
and precessing spin distributions has important consequences for spin dynamics and spin-related 
effects in general. In the latter part of the article extrinsic effects such as skew scattering and side 
jump will be discussed, and it will be shown that these effects are also modified considerably by 
spin precession. Theoretical and experimental progress in all areas will be reviewed. 
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^ I. INTRODUCTION 

Spin electronics seeks to harness the spin degree of freedom of the electron, aiming to generate and maintain a 
spin polarization for possible use in information storage through the manufacture of magnetic memory devices. For 
spin-based quantum computing the manipulation of a spin polarization is also crucial. A spin polarization can be 
genegrated by optical, magnetic and electrical means. Optically, one uses light to transfer angular momentum to 
charge carriers during interband transitions. [1| This reliable procedure has been tested successfully decades ago, yet 
devices that require an input of light are currently impractical. Magnetically, one uses the Zeeman effect, to which 
similar observations apply. Electrically, one drives a current from a spin-polarized material into an unpolarized one and 
thus carry spin-polarized electrons across the interface. This works in the case of metals but modern technology relies 
£j | overwhelmingly on semiconductors, and spin injection from ferromagnetic metals into semiconductors is hampered by 
the resistivity mismatch Q at the interface between the metal and the semiconductor, which causes most of the spin 
polarization to be lost at the interface. Another avenue explores the use of ferromagnetic semiconductors [3] as spin 
injection devices, but room temperature ferromagnetism in semiconductors remains a long-term goal. 
■ As a result of these challenges, purely electrical means of spin generation have been sought. The past decade 
7-H , has witnessed an explosion of interest in electrically-induced spin phenomena in semiconductors and metals, which 
7-H ■ has accompanied significant experimental and theoretical progress. It has been found that the application of an 
electric field gives rise to a nonequilibrium spin density in the bulk of the sample, discussed in Refs. and a 

\Q • nonequilibrium spin current, discussed in Refs. The appearance of a steady-state spin density in an electric 

field was predicted several decades ago [|| and subsequentl y o bserved in tellurium. [5[ This work was followed by a 
0^ , number of theoretical @, H, and experimental [ToL [lTlTl2| studies. Spin currents also have a long history, their 
theoretical prediction stretching back almost four decades to the work of Dyakonov.[l3j The theoretical work was, 
, however, not followed by experiments, and the topic was silent for over a quarter of a century, until the theory picked 
up again, led by Hirsch [l4 | , Murakami et al. [HI and Sinova et al. fl6j During the recent surge of interest in electrically- 
induced spin phenomena the focus has been principally on spin currents, with the ultimate goal of measuring spin 
flow in a given direction. In the absence of a spin-ammeter, initially the only way of detecting a spin current relied on 
measuring the spin accumulation it produced at the edge of the sample. This procedure, as will be discussed below, 
is complicated by the fact that the spin current is often not well defined and its relationship to spin accumulation 
is not clear. Thus one can detect a spin accumulation at the edge of the sample as a result of the spin current, but 
inferring the size of the spin current from such a measurement is beyond current understanding. Nevertheless, in just 
a few years experiment has been surging ahead. First, in a recent experiment, Crooker and Smith [l7| imaged the 
flow of spin in semiconductors in the presence of electric and magnetic fields, as well as strain. Following that, spin 
accumulations resulting from spin currents were measured by several groups. [H, [H, [2(| HH, [H, [23], [24[ In what is 
surely a new generation of experiments, Valenzuela and Tinkham[25j made use of the fact that a spin current gives 
rise to a transverse charge current. Relying on a nonlocal technique, their group successfully detected a charge current 
flowing as a result of a spin current in Al metal. Cui et al. [26] also observed an electrical current induced by a spin 
current, though in that case the spin current was generated optically. This technique was used by several groups 
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[13, [H, [H, [33, [M| shortly afterwards in the study of platinum- and gold-based structures. In gold an unusually large 
spin current was observed[3(| even at room temperature. Room temperature spin currents were also reported by 
Stern et al 21] in ZnSe. Thus, from their initial observation four years ago, spin currents have begun an experimental 
revolution, and the field has matured to the point where experiment and theory interact constructively. 

The phenomena discussed in this review are all related to spin-orbit interactions, which are present in the band 
structure and in potentials due to impurity distributions. Band structure spin-orbit interactions are frequently the 
most important factor determining spin dynamics in solids, and are particularly important in semiconductors, where 
carriers involved in steady-state processes have wave vectors well within the first Brillouin zone. At these wave vectors, 
while the spin-orbit contribution to the carrier energy is still smaller than the kinetic energy, it plays an important role 
in determining the energy spectrum. Band structure spin-orbit coupling may arise from the inversion asymmetry of 
the underlying crystal lattice (86j , from the inversion asymmetry of the confining potential in two dimensions [87| , and 
may be present also in inversion symmetric systems. [88| Spin-orbit coupling is also in principle present in impurity 
potentials. Spin-orbit interactions in the impurity potentials cause skew scattering, or asymmetric scattering of up and 
down spins, giving scattering-dependent or extrinsic contributions to spin currents. Finallly, spin-orbit interactions 
cause a modification of the position and velocity operators, which affects the interaction with an electric field as well as 
the scattering term. This last mechanism is referred to as side jump and is typically classified as extrinsic although its 
contributions are manifold. These extrinsic mechanisms were discovered in the anomalous Hall effect [89l [90I . Wu , |92| 
but have acquired a new relevance in recent years in electrically-induced spin phenomena. 

In charge transport, the steady state is characterized by a nonequilibrium density matrix that is divergent in the 
clean limit, indicating a competition between the electric field, accelerating charge carriers, and scattering, which 
inhibits their forward motion. The density matrix is a scalar and the correction to it caused by the electric field is 
cx , where rii is the impurity density. When considering the steady state of systems with spin-orbit interactions 
three different situations are distinguished: the case when only extrinsic mechanisms are important, the case when 
only band-structure mechanisms are important, and the case when both extrinsic and band-structure mechanisms are 
important. The differences between these cases are significant and nontrivial. 

If band structure spin-orbit interactions are negligible and only extrinsic mechanisms are important, there is no 
spin precession and the correction to the spin density matrix due to an electric field is very similar to the correction to 
the scalar density matrix in charge transport. This correction is oc n~ x and it gives rise to a steady-state spin current 
but not to a steady-state spin density. There have been many studies of extrinsic spin currents in semiconductors 
(Refs. [[Z3-[7|]) and metals (Refs. 

Nonequilibrium corrections that arise as a result of band structure spin-orbit coupling in crystal Hamiltonians 
represent a different kind of interplay between the electric field and scattering processes. The spin-orbit splitting 
of the bands gives rise to spin-dependent scattering even from spin-independent potentials, and spin currents and 
spin densities in an electric field arise from linearly independent contributions to the density matrix. The steady- 
state density matrix contains a contribution due to precessing spins and one due to conserved spins. Steady-state 
corrections oc n~ x are associated with the absence of spin precession and steady-state corrections independent of 
are associated with spin precession. Steady-state corrections oc n~ give rise to spin densities in external fields while 
steady-state corrections independent of m give rise to spin currents in external fields. Scattering between these two 
distributions induces significant corrections to steady-state spin currents, and may cause spin currents to vanish. 

The phrase spin current refers to the flow of spins across a sample and, if spin were conserved, one could distinguish 
between spin-up and spin-down c harg e currents. In systems with band structure spin-orbit interactions the spin 
current is not well defined [32l . l33l . l34j and its relationship to spin accumulation is unclear. Spin transport in these 
systems usually does not involve charge transport as the charge currents in the direction of spin flow cancel out. 
In addition, since spin currents may be accompanied by steady-state spin densities, this makes it more difficult to 
separate experimentally the signals due to these contributions. The study of spin currents in systems with band 
structure spin-orbit interactions has thus encountered a number of profound physical issues on which no consensus 
exists at present. There has been considerable research has on definitions of spin currents in such systems (Refs. [f32|- 
[35|]). the role played by the symmetry of the underlying crystal lattice. |3q] whether spin currents as a result of 
band structure spin-orbit coupling are transport or background currents, [371 ] the relationship between band structure 
spin currents and spin accumulation, (Refs. [[Hj]-[4i|]) and the form of the Maxwell equations in systems with band 
structure spin-orbit coupling. (4(| H?| Most theoretical studies have focused on metals, common semiconductors and 
asymmetric quantum wells with band structure spin-orbit interactions, such as Refs. [ |16| - [76j|] . 

When both band structure and extrinsic mechanisms are present, the situation becomes considerably more compli- 
cated. There is no simple association between spin conservation and spin densities, or between spin precession and 
spin currents. Both spin densities and currents may contain terms oc n~ l and independent of m. A careful analysis 
shows that the presence of band structure spin-orbit interactions causes skew scattering and side jump to behave very 
differently in the steady state, and under certain circumstances, the contributions of these extrinsic mechanisms to the 
spin current vanish. The interplay of band structure and extrinsic mechanisms has been considered in Refs. [[80j-(83|]. 
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This review will cover steady-state spin densities and currents within the framework of a density matrix theory. 
Such a framework is important in order to demonstrate the unity behind all observed phenomena. I will consider 
large, uniform systems, and work in momentum space. Although many observations in this entry are general, the 
discussion will focus on non-interacting spin-1/2 electron systems, which are pedagogically easier. This review is based 
on a kinetic-equation formalism equivalent to linear response theory and correspondences with other methods will be 
identified, such as linear response theories based on Green's functions and semiclassical wave packet dynamics. 

This review article is structured as follows. Section 2 will introduce a density matrix picture of spin dynamics, briefly 
sketching the way different theories stem from the quantum Liouville equation. In section 3 a brief introduction to 
the spin-orbit interaction in spin-1/2 systems is given, after which spin densities and spin currents as a result of band 
structure spin-orbit interactions are discussed, as well as their relationship to spin precession. Section 4 is devoted 
to spin densities and spin currents that arise as a result of extrinsic mechanisms, together with the case when both 
extrinsic and band-structure mechanisms are present. Section 5 discusses briefly various definitions of spin currents 
and the relationship between spin currents and spin acumulation. Section 6 outlines the role of the underlying crystal 
lattice in the establishment of spin densities and currents and section 5 is concerned with open issues such as the 
definition and nature of spin currents, as well as the relationship between spin currents and spin accumulation. In 
section 7 the experimental situation is summarized and in closing future directions are discussed. Related topics such 
as the quantum spin-Hall effect [93| and the spin-Hall insulator [94[ are beyond the scope of this review. 

II. DENSITY MATRIX PICTURE OF SPIN DYNAMICS 

Electrically-induced spin phenomena encompass a wide variety of processes that have been studied using different 
methods. In order to bring out the unity behind the processes involved, as well as behind the theoretical approaches 
employed, it will be useful to have a unified framework in which electrically-induced spin phenomena can be discussed. 
In this section such a framework will be constructed beginning with the quantum Liouville equation. This equation, 
in one form or another, is the starting point of most theories of spin dynamics, and this section will outline the way 
various approaches are related to each other. The focus will be on systems with large homogeneous systems with long 
mean free paths, and diffusion will not be considered explicitly. Electric fields will be assumed uniform. 

A. Quantum Liouville equation 

A system of non-interacting spin-1/2 electrons is represented by a one-particle density operator p. The expectation 
value of an observable represented by a Hermitian operator O is given by Tr(/SO), where Tr denotes the most general 
operator trace, which involves summation over discrete degrees of freedom and integration over continuous ones. The 
usual matrix trace will be denoted by tr. The dynamics of p are described by the quantum Liouville equation, 

^ + ~ [H + U dls + eE ■ f, p] = 0. (1) 
at h 

The Hamiltonian H contains contributions due to the kinetic energy and spin-orbit coupling, while r is the position 
operator. The effect of the lattice-periodic potential of the ions is taken into account through a replacement of 
the carrier mass by the effective mass. The potential U dls accounts for scattering processes, which may be due to 
impurities, phonons, surface roughness, or other perturbations. This review focuses on impurity scattering, as the 
effects discussed are frequently observed at low temperatures, where scattering due to phonons may be neglected. It 
is assumed that the electric field is small and a solution is sought to first order in the electric field. The solution of the 
Liouville equation ([T]) in the absence of the external field is provided by pa , the Fermi-Dirac function. The correction 
due to the external field, pe, satisfies 

*+ J i*+e-, w --«(*.»,«. ( 2) 

Equation ([2]) is projected onto a complete set of time- independent states of definite wave vector {\ks)}. These states are 
not assumed to be eigenstates of H. The matrix elements of p in this basis will be written as pkk' = Pkk' = (ks|/5|fcV), 
with corresponding notations for the matrix elements of H and U dls . Spin indices will not be shown explicitly in the 
subsequent discussion, and it will be understood that the quantities pkk' , Hkk' , and Uj^, are matrices in spin space. 
Pkk 1 is referred to as the density matrix. With our choice of basis functions of definite wave vector, matrix elements of 
the Hamiltonian H^k' = Hk &kk' are diagonal in k. However, since the Hamiltonian contains spin-orbit coupling terms, 
matrix elements Hk are generally off-diagonal in spin space. Matrix elements of the scattering potential U^, are 
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off-diagonal in k. Matrix elements diagonal in k in the scattering potential would lead to a redefinition of Hk, which 
is analogous, in Green's function formalisms, to the offset introduced by the real part of the self energy. Scattering 
is assumed elastic and impurities uncorrelated, and the normalization is such that the configurational average of 
(ks\U \k's')(k's'\U \ks) is (rii\Ukk'\ 2 Sss')/V, where nt is the impurity density, V the crystal volume and Ukk 1 the 
matrix element of the potential of a single impurity. Configurational averages over terms of higher order in jj dls are 
performed in a similar fashion. [96[ It is assumed that epr/h^ 1, where Ef is the Fermi energy and r a characteristic 
scattering time. This is equivalent to the assumption that the carrier mean free path is much larger that the de 
Broglie wavelength. None of the methods presented below are valid once epTp/h becomes comparable to unity. 

Spin densities and spin currents are obtained as expectation values of spin and spin current operators. The spin 
operator is given by s CT = (h/2) <r a , where a a is a Pauli spin matrix. The spin current operator will be taken to be 
J? = (1/2) {s CT , v 1 }, where the velocity operator is v l — (1/fr) dHk/dki. Both operators are diagonal in k. One is 
therefore primarily interested in the part of the density matrix diagonal in k. With this observation in mind, pkw is 
divided into a part diagonal in k and a part off-diagonal in k as pkw = fk Skk' + 9kk' , where, in it is understood 
that k =/= k' . Further, fk is decomposed into a scalar part and a spin-dependent part, fk — «fc 1 + Sk, with 1 the 
identity matrix. In an electric field, fk = fok + $Ek, where fEk is a correction linear in E and has a corresponding 
decomposition into a scalar part and a spin-dependent part. At the end one takes the trace of the spin and spin 
current operators with fEk given by Eq. Q. The final result is usually expressed in many ways. 

It should be noted that the electric field induces coherence between bands. In equilibrium the Fermi-Dirac function 
fok commutes with the Hamiltonian and is stationary. However the source term due to E in Eq. @ does not in 
general commute with the Hamiltonian. If the problem is considered in the basis of eigenstates of Hk, the source 
term couples different energy bands. If a basis is used in which one spin component is a good quantum number, then 
the source term couples up and down spins. 



B. Linear response Green's function formalism 

The linear response Green's function formalism is closely related to the density matrix formalism but relies on a 
somewhat different way of regarding the problem [961 ]. The electric field-induced correction to the density matrix pe 
is found by applying the time evolution operator corresponding to the total Hamiltonian H + U dls , incuding disorder, 
in Eq. ([2]). The matrix elements of the time evolution operator in the basis {|fcs)} constitute the Green's function 
G s k i(t), defined by 

Gii(t) = (ks\e-^+° d,s ^\k's'). (3) 

Spin indices are again suppressed. All the information about the system is contained in the Green's function, which 
is also referred to as the propagator, since it gives the amplitude for a particle to propagate in time t from state 
\ks) to state |fcV). The Green's function is also the kernel of the Schrodingcr equation and represents the response 
of the system to a perturbation that is localized in real space or momentum space. The total perturbation is built 
up as a sum of localized ones. The problem can be formulated in any basis but for our purposes it is still the most 
convenient to work in the basis of eigenstates of H so k- In general the Green's function is a matrix in spin space. For 
practical purposes one also defines the retarded and advanced Green's functions, G R and G A , as G R — —i G9(t) and 
G A = iG6(—t), with 8 the Heaviside step function. The electric field-induced correction to the part of the density 
matrix diagonal in k is 

f Ek = - / dt'G* k ,(t')[r,p } k ,G A , k (-t') (4) 

n k> 

At this stage it is easiest to carry out a Fourier transformation with respect to time. The integration in Eq. ^ has 
the same form except the time variable is replaced by the frequency variable u> 

fEk = --^--Y. / duG* k ,{u)[r,p } k 'G£, k {-u). (5) 
n k< J -°° 

The correction fEk to the density matrix needs to be averaged over impurities, and it is assumed henceforth that 
only the impurity averaged density matrix is of interest. We denote this impurity average by the symbol (), but the 
impurity-averaged density matrix will be abbreviated as(/Efc) = fEk- The full Green's function is determined by the 
total Hamiltonian H + U dls including the scattering potential. The Green's function is expanded in the scattering 
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potential, and one is only interested in the Green's function averaged over impurity configurations. For uncorrelated 
impurities the impurity-averaged Green's function, denoted by (G), obeys the recursion relation 

{G R / A {k,u)) = G^ /A {k,Lu) + G^ /A (k,u;)^ A (k,Lu)G^ A (k,Lu) (6) 

in which Go is the equilibrium Green's function, corresponding to the Hamiltonian H without disorder, and the self 
energy £ is typically evaluated in the first Born approximation 

s«M(fc, w ) = mj0yi \U kk ,\ 2 G«/ A {k',u) (7) 

with d the dimensionality of the system. The correction f E k to the density matrix depends on the product of two 
Green's functions. After the average over impurity configurations we obtain terms which are expressible as a product 
of two individual impurity-averaged Green's functions and cross terms, which connect different Green's functions, and 
which are not expressible as a product. These form a vertex Tkk 1 defined by 

(G R (k, oj)G A {k', u)) = (G R (k, W )) (G A (k, w)) + T kk , (to) {G R {k', co))(G A (k\ w )>. (8) 

The vertex is in general expressed as a sum of different classes of diagrams, of which the most relevant to transport 
are the ladder diagrams and the maximally crossed diagrams. It also satisfies the recursion relation 

T kk ,(io) = M kk ,{uj) + J ^ 1 M kk , l (Lu)(G R (k'',u))(G A (k' , ,u)T k ,, k ,(Lu)) (9) 

where M kk i [to) can also be expressed in terms of diagrams as shown in Ch. 8 of Ref. [[96|]- Once the impurity-averaged 
Green's function (G) and the vertex are known, f E k can be found. It will be noticed that in the linear response Green's 
function approach disorder is contained in the self energy E(fe) and in the vertex T kk i . The response function contains 
all the disorder up to the desired approximation. Once the Green's function and the vertex function are found, the 
trace with the spin and spin current operators is taken. 



C. Kinetic equation formalism 



The kinetic equation formalism can be derived from the quantum Liouville equation, as will be done below, or using 
a Keldysh Green's function formalism. [EtJ The kinetic equation obtained is naturally the same. Starting from the 
quantum Liouville equation, this can be broken down into equations for f k and gkk' 



dl 



Ek 



dt 

dgEkk' 



[H,g E . 



Ek\ 



kk' 



T [U, g E ]kk - - [eE ■ r, p Q } k k, 



[U,fE+9E]kk> 



dt h l '^ J h 

In the first Born approximation the solution to Eq. (jlObp can be written as 



gEkk' 



dt , e -Wt'/H 



u,f E (t-t') 



(10a) 
(10b) 

(11) 



Since EFT p /h 3> 1, we shall expand f(t — t') in the time integral around t and, noting that terms beyond f(t) are of 
higher order in the scattering potential, we shall only retain the first term, fit). The equation for f k then becomes 



~dT h [Hfc ' fEk] + J{lEk) = J^\-dk- ^' /ofc] 



in which the scattering term J(f k ) is given by 



J [fEk) - / 

ft Jo 



dt' 



ft -iHt'/H 



U,f E (t) 



jiit'/h 



kk 



(12a) 



(12b) 



and the source term on the RHS contains the covariant derivative with respect to fc, which takes into account the 
fact that the basis functions themselves may depend on k. It includes the gauge connection matrix 1Z = TLss' = 



6 



(ks\id/dk\ks'). The response function, as well as the scattering term need to be averaged over impurity configurations, 
and the notation f^k is used henceforth as above to denote the impurity-averaged density matrix. The scattering 
term can be expanded further in the basis in spin space spanned by spin eigenstates | |) and | J.) (the Pauli basis.) 
In this basis the connection matrix 7L vanishes and the scattering potential Ukk' = Ukk' 1 is diagonal in spin space. 
Converting sums over wave vector into integrals ^2 k , — ► V J ^ i an< i performing the time integral the scattering 
term can be expressed in the form (J + J&) (n k ) + Ja{Sk), with 

2-7m f d d k' 

Ja(Sk) = J j^y[\Ukk'\ 2 {Sk- S k ')S(e ok - £ok>) (13a) 

- . . l-KUi f d d k' 2 1 . d ... 

Jb(nk) = — r— / 7^—rj\Ukk>\ (nfe - nfeO - cr • (S2fe — Slfe')- d{e k ~ eofc'J- 

h J (2n) d 2 cteofc 

The term Jf, (rife) in Eq. (|13bj) illustrates the fact that, when spin-orbit interactions are present in the band structure, 
even spin-independent scattering potentials give rise to spin-dependent terms in the scattering integral. For potentials 
diagonal in spin space and spin-degenerate bands, Eq. (|12bjl simplifies to the customary Fermi's golden rule. Therefore, 
Eq. (|12p is a generalization of Fermi's golden rule that explicitly takes into account the spin degree of freedom. 



D. Boltzmann-wave packet formalism 



The density matrix is the most complete description of a physical system. It contains all the relevant physics, 
including interband coherence due to the electric field and scattering. It represents the system as a whole and 
accounts for the individual particle motion, which occurs between collisions, as well as for the distribution of electrons 
in phase space and the way it is altered by collisions. As was shown above, linear response theories based on Green's 
functions and theories based on kinetic equations originate directly from the Liouville equation for the density matrix. 
A third approach, which may be called a semiclassical or Boltzmann-wave packet approach, separates the motion 
of the charge carriers from their distribution in phase space. To determine the dynamics of the charge carriers one 
envisages them as being represented by wave packets, and the phase space distribution of wave packets is given by a 
function of the Boltzmann form. One use of the word semiclassical, which will be adopted in this work, is to refer to 
theories which consider the position and momentum of a particle simultaneously. Semiclassical approaches exploit the 
smooth variation of transport fields on atomic length scales in order to provide intuitive descriptions of steady-state 
processes, which can be easily extended to cover inhomogeneous systems and spatially dependent fields. 

The semiclassical approach is also closely related to the density matrix and the Liouville equation. Firstly, in order 
to separate the particle dynamics and the phase space distribution one must be able to label the particles by a band 
index, therefore the theory must be formulated in the basis of eigenstates of H k . Since these eigenstates are usually 
fc-dependent it is the covariant derivative that enters the source term in the kinetic equation. The point at which 
semiclassical theory branches out of the kinetic equation is Eq. (|12ap . In order to turn the density matrix into a 
Boltzmann distribution function one takes the diagonal elements of this equation. If we consider for simplicity a spin- 
1 /2 system, which is made up of two bands that will be called 1 and 2, we can label the diagonal elements of the density 
matrix as fik and fek- These become the Boltzmann distributions for bands 1 and 2. The commutator [Hk, fk] has 
no diagonal elements and J(fk) reduces to one scattering term for each band, which is given by Fermi's Golden Rule, 
plus interband scattering terms which couple fik and f2k- In this way Eq. (|12a[) reduces to two coupled Boltzmann 
equations for the two bands and the phase space distribution has been separated out of the kinetic equation. 

In reducing the kinetic equation to a series of Boltzmann equations one neglects interband coherence due to the 
electric field, and this must be recovered elsewhere. In the semiclassical approach, this occurs in several ways. Firstly, 
one needs to consider the nature of the carriers in band s and invoke explicitly the fact that carriers are described 
by wave packets. This is the point at which the formalism turns into a true semiclassical theory. The construction 
of a wave-packet representing a charge and spin carrier in band s, which has real and /c-space coordinates (r cs , k cs ), 
has been thoroughly treated by Sundaram and Niu [97| and its extension to multiple bands was carried out in Ref. 
[[98j], The macroscopic distribution associated with a quantum mechanical operator O is no longer given simply by 
the expectation value of the operator O, but involves explicitly a sum over wave packets centered at (t* cs , k cs ). For 
example the spin density distribution is 

5 CT (r,i)=^ J d 3 k cs J d 3 r cs f s (k cs ,t)(S(r-r)s°) s , (14) 

where the bracket indicates quantum mechanical average over the wavepacket with charge centroid (r cs ,k cs ). An 
analogous expression exists for the spin current distribution. The wave-packet expectation values () s in Eq. (|14|) 
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FIG. 1: For a particle of finite extent the charge and spin distributions in real space are in general do not coincide. The same 
is true of the charge and spin distributions in reciprocal space. 



eventually yield functions of (r cs ,k cs ), the dynamics of which are discussed below. It may appear contradictory that 
one has to integrate over the real-space coordinates of the wavepackets r cs even in the case of homogeneous systems 
studied in this work, where the distribution function is not a function of r cs . This is because even in a homogeneous 
system of the center of charge and the center of spin are not the same, so variations of the spin distribution on the 
spatial scales comparable to that of the wave packet need to be taken into account, as shown in Fig. 1. All these 
terms must be considered in order to reach agreement with approaches based directly on the Liouville equation. 

In a constant uniform electric field E the coordinates of the wave packet center (r cs ,k cs ) drift according to the 
semiclassical equations of motion [97j 



-eE 



de. 



foes = heExf s 

<9k c 



(15) 



where J- a = xTZ ss represents the Berry, or geometrical curvature [971 ] . with the gradient operator in momentum 
space. A careful analysis reveals that the Berry curvature also represents coherence between bands and must be taken 
into account in order to obtain agreement with approaches based directly on the Liouville equation. It is emphasized 
that the Berry curvature is not a result of the particles being described by wave packets, which have a finite extent in 
real and momentum space. Rather, it emerges from the phase of the Bloch wave functions involved in the construction 
of the wave packets, implying that this formulation of semiclassical transport theory is a useful tool for capturing 
Berry-phase effects. 

The electric field and disorder potential also give rise to a nonadiabatic mixing of the bands such that the basis 
states \ks) become \ks), given by 



\ks) = \ks) 



E 

k' ,s'^s 



[eES, 



&ks &ks 



\ka'), 



(16) 



where the \ks) are the unperturbed eigenstates. The \ks) also form a complete set. The distribution functions 
f s (k ca , t) are made up of an equilibrium part and a part linear in the electric field, f s (k cs , t) — fa s (k cs , t) + /E S (fe cs , i), 
and the expectation values () s in Eq. (|14[) also have terms of zeroth and linear orders in the electric field. Therefore 
to first order in the field the expectation value of operator O contains terms arising from the product of fEs(k cs ,t) 
with the zeroth order term in () s as well as from the equilibrium fo s (k cs , t) multiplied by the term in () s linear in the 
electric field. 

The above exposition shows that interband coherence lost in passing from the density matrix to the Boltzmann 
approach is thus recovered in many ways. It is present in the Berry curvature, in the nonadiabatic correction to 
the wave functions, and in the interband scattering terms coupling the two distribution functions. The construction 
of a correct semiclassical theory that contains all the physics relevant to transport is seen to entail a sizable effort, 
particularly in the case of transport of non-conserved quantities such as spin. 
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III. BAND-STRUCTURE SPIN DENSITIES AND CURRENTS 



This work will now focus on the kinetic equation formalism and determine the spin densities and spin currents 
induced by external electric fields under a variety of circumstances. It is most straightforward to begin with effects 
originating in the band structure. The Hamiltonian of spin-1/2 electron systems typically contains a kinetic energy 



term and a spin-orbit coupling term, H k = -§^V + Hff, where m* is the electron effective mass. In spin-1/2 electron 
systems, band structure spin-orbit coupling can always be represented as a Zeeman-like interaction of the spin with 
a wave vector-dependent effective magnetic field fl k , thus Hff — {ft/2) cr ■ $l k . Common examples of effective fields 
are the Rashba spin-orbit interaction, 18 7| which is often dominant in quantum wells with inversion asymmetry, and 
the Dresselhaus spin-orbit interaction, [86f which is due to bulk inversion asymmetry. 

An electron spin at wave vector k precesses about the effective field £l k with frequency £l k /h = \Sl k \/h and is 
scattered to a different wave vector within a characteristic momentum scattering time r. Within the range eft/H^> 1, 
the relative magnitude of the spin precession frequency £l k and inverse scattering time 1/t define three qualitatively 
different regimes. In the ballistic or clean regime no scattering occurs and the temperature tends to absolute zero, 
so that Eft/H — > oo and fl k r/h — * oo. The weak scattering regime is characterized by fast spin precession and little 
momentum scattering due to, e.g., a slight increase in temperature, yielding eft/H 3> fl k r/h 3> 1. In the strong 
momentum scattering regime Eft /ft > 1 > fife t/H. 

The kinetic equation (|12ap has a scalar part which determines riEk 

dnEk , * , . eE dn ok 

-g r + Jo(n Ek ) = — -^. (17) 

The solution of this equation is given by the well-known expression 

eEr p dn ok 

nEk = — - "air' (18) 

in other words, riE k describes the shift of the Fermi sphere in the presence of the electric field E. The expression 
for the momentum relaxation time t p is a little different depending on the dimensionality of the system. Using 7 to 
denote the relative angle between k and k' , 



1 riimk 
t p 2Trh 3 

riim 



[ d~f\U kk ,\ 2 sin 7(1 -COS7) (3D) (19a) 
Jo 

/^7|£4 fe '| 2 (l-cos 7 ) {2D). (19b) 



The spin-dependent part of the nonequilibrium correction to the density matrix SEk is the spin density induced by 
E. Its time evolution is governed by SE k is 

—q7- + t l H k, S Ek \ + J {S Ek ) = -z 7^7- ~ Jb [n Ek )- (20) 

Spin-dependent scattering gives rise to a renormalization of the driving term in the equation for SE k with no analog 
in charge transport, see Eq. (fTT)) . 

Since an electron spin at wave vector k precesses about Q k , the spin can be resolved into components parallel and 
perpendicular to £l k . In the course of spin precession the component of the spin parallel to £l k is conserved, while the 
perpendicular component is continually changing. It will prove useful in our analysis to divide the spin distribution 
into a part representing conserved spin and a part representing processing spin. The effective source term, which 
enters the RHS of Eq. (j20|) . is divided into two parts, {eE/h) ■ dSo k /dk — Jb {riE k ) = ^Ek\\ + ^Ek±- ^Bk\\ commutes 
with Ht° while T.E k ± is orthogonal to it and tv{Y^E k ± Hf?) — 0. Projections onto and orthogonal to are most 
easily carried out by defining projectors P|| and P± by their actions on the basis matrices Oi as described in Ref. [[36]]. 
SEk is likewise divided into two terms: Sem, commuting with the spin-orbit Hamiltonian and SE k ±, orthogonal to 
it. S E k\\ is the distribution of conserved spins while SEk± is the distribution of precessing spins. Equation (|20[) is 
divided into separate equations for SE k \\ and Ssk± 

" ! ' L ' V <P ll Jo{S E k) = ZEk\\, (21a) 



at 

dS Ek ± i_ 
dt ft 



[H k ,S Ek ±] + P±Jo {S Ek ) = ^ Ek ±- (21b) 
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The absence of the commutator [Hk, Sem] — in Eq. (|21a|) indicates the absence of spin precession, while the 
commutator [Hk, SEk±] in Eq. (|21b[) represents spin precession. In order to solve Eqs. (|21a[) and (|21b[) for arbitrary 

scattering, it is necessary to expand SEk\\ and SEk± in "i, as SEk\\ = S EI ^ + S E ^ + 0(n,i) and SEk± = S E l ± + 0(ni). 
This is an expansion in the parameter h/(Qk T p) and is most suited to systems in the weak scattering regime. The 
expansion of SEk\\ starts at order —1, a fact which can be understood by inspecting Eq. (|21ap . In the steady state 

the time derivative drops out, and the operator Jo is first order in rii, while the right-hand side is independent of rn. 
As a result, the expansion of the solution must start at order —1. Equation (I21b[) for SEk± tells us that, since Hk is 
independent of rii and the right hand side is also independent of rn, the expansion of S^fex must start at order zero. 
Equation (|21a[) can be solved iteratively for any scattering 

S£J = £ £fc ||T + P|, (j^ J d9'\U k k> I 2 £ £fc ||)r 2 + • ■ • (22) 

where to = riim* J ' d9'\Ukk'\ 2 /(^ttH 3 ) is the quantum lifetime of the carriers. Above and henceforth integrals over 
wave vectors will be represented as two-dimensional, and 9' will refer to the polar angle of k' . The extension to three 
dimensions is straightforward. The equations for higher orders in rii are easily deduced. However, the term of order 
— 1 is by far the dominant one in the weak momentum scattering regime and is expected to be dominant over a wide 
range of strengths of the scattering potential. 

It is evident that the steady state for conserved spins involves no spin precession, and that the correction Sem 
depends explicitly on the nonequilibrium shift in the Fermi surface. In addition, scattering terms contain only the 
even function \lAkk 1 \ 2 ■ As a result, Sem does not give rise to a spin current. Inspection of Eq. (|22[) shows that integrals 
of the form J dd J? Sem contain an odd number of powers of k and are therefore zero. In the absence of impurity 
spin-orbit interactions, the distribution of conserved spins can give no spin current. It can, however, give rise to a 
nonequilibrium spin density since integrals of the form J d9 s a SEk\\ contain an even number of powers of k and may 
be nonzero. 

The leading term S E l ± is found to be 



c (o) _ ft t ■ O fc x [S E k± - P±Jp (SEk\\,> v> 
^ fc ^ _ 2 W k ' 1 > 

where we have written T*Ek± = (1/2) Sejji • a and Sem = (1/2) Sem ' <T - The result expressed by Eq. (|23p is valid 
for any elastic scattering. An argument similar to that given above shows that SjSl ± cannot lead to a nonequilibrium 
spin density (although, as will be shown below, higher-order terms in SEk± can contribute to the spin density). For, 
taking the expectation value of the spin operator, one arrives at integrals of the form J dO s a S E lj_, which involve odd 
numbers of powers of k and are therefore zero. This term in the distribution of precessing spin does, however, give 
rise to nonzero spin currents, since integrals if the form J dO J? SEk± contain an even numbers of powers of k and 
may be nonzero. Consequently, in the absence of spin-orbit coupling in the scattering potential, nonequilibrium spin 
currents arise from spin precession. 

The dominant contribution to the nonequilibrium spin density in an electric field exists because in the course of spin 
precession a component of each individual spin is preserved. For an electron with wave vector fc, this spin component 
is parallel to f2fc. In equilibrium the average of these conserved components is zero. However, when an electric field is 
applied, the Fermi surface is shifted, and the average of the conserved spin components may be nonzero, as illustrated 
in Fig. 2. This intuitive physical argument explains why the nonequilibrium spin density oc r" 1 and requires scattering 
to balance the drift of the Fermi surface. It is interesting to note, also, that, although spin densities in electric fields 
require the presence of band structure spin-orbit interactions and therefore spin precession, the dominant contribution 
arises as a result of the absence of spin precession. Band structure spin currents on the other hand are associated 
with displacement of spins. The relation between spin currents and spin precession was made explicit by Sinova et al. 
[l6| Both Sem an( i SEk± are invariant under time-reversal. As a result, the tensor characterizing the response of 
spin currents to electric fields is invariant under time reversal, whereas the tensor characterizing the response of spin 
densities to electric fields changes sign under time reversal. 

It will be noticed that the driving term in the equation for the nonequilibrium spin distribution SEk is renormalized 
by the term J s (riEk), which accounts for spin-dependent scattering. In addition, Eq. (|21b[) shows that scattering 
mixes the distributions of conserved and precessing spins. When one spin at wave vector k and precessing about 
Jlfe is scattered to wave vector k' and precesses about its conserved component changes, a process which 

alters the distributions of conserved and precessing spin. Consequently, scattering processes in systems with spin- 
orbit interactions cause a renormalization of the driving term for the spin distribution as well as scattering between 
the conserved and precessing spin distributions. Contributions due to these two processes are contained in vertex 
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FIG. 2: Effective fieid ftk at the Fermi energy in the Rashba modei (a) without [E = 0) and (b) with (E > 0) an external 
electric field. 



corrections to spin-dependent quantities found in Green's functions formalisms. For short-range impurities it is easy 
to show [36l.r75| that the scattering correction P±Jq {SEk\\) to Eq. (|23p depends on the steady state spin density, and 
it becomes evident that the existence of a nonzero nonequilibrium spin density does affect the spin current. This fact 
was first pointed out by Raimondi [75[ for the Rashba model. 

The spin-Hall current was initially determined for spin-3/2 holes in GaAs [la ] and for an asymmetric quantum well 
[l6| in which the spin-orbit interaction is described by the Rashba model. The latter calculation yielded a spin-Hall 
conductivity a% y — e/(87r), in which <7* fe is understood as referring to spin component i flowing in direction j in 
response to an electric field applied along k. However, it was subsequently shown that a more careful treatment of 
disorder renormalizes this result to zero. In fact in two dimensions, for Hamiltonians linear in wave vector, the spin 
current vanishes for short-range impurities (Refs. [ 50]- 55]]), for small-angle scattering [56H57| and in general for any 
elastic scattering. [33] However, it does not vanish in a generalized Rashba model as was shown by Krotkov and Das 
Sarma. |68l| For spin-orbit Hamiltonians characterized solely by one angular Fourier component N the spin current 

octv. da] 

Ab initio calculation of band structure spin currents were also performed by Guo et al. [72J in semiconductors and 
by Yao et al.ff^ in semiconductors and simple metals such as tungsten, platinum and gold. Spin transport in metals 
usually requires complex band structure calculations, which are typically done numerically. Yao et aL [731] found that 
the band structure spin-Hall effect in metals can be significantly larger than in semiconductors and that the spin-Hall 
conductivity can even undergo sign changes under certain circumstances. 



A. Spin adiabatically following a magnetic field 

The existence of a spin current in an electric field and its association with spin precession, can be understood from 
a straightforward argument due to Sinova. [T(| A spin precesses about an effective magnetic field f2fc which depends 
on k. The electric field changes the wave vector k and in the process changes H^, so the spin is now precessing 
about magnetic field that is changing adiabatically. In the adiabatic limit the spin follows the magnetic field. [§!| In 
particular, in 2D if the magnetic field is in the plane of the spin, the spin never acquires a significant out-of plane 
component. However the spin in general acquires a small out-of-plane component. Consider a generalized magnetic 
field O and take the simple example in which the magnetic field || x and has a small y-component which increases 
linearly with time, such that £l y — et. The spin starts out along x, initially parallel to the field. In both cases, s x will 
be considered large compared to s y , s z and in this section h — 1. The adiabatic condition means that the magnitude 
of the magnetic field changes little over one revolution of the electron spin about it. This means that, if we start 
at t = and consider a small time increment t, f2(r) — O(0) << O(0) . 

O(r)=O(0) + ^| T 

(24) 

O(r)-O(0) = ^| r^^| r«O(0). 

The adiabatic limit in this case is the limit of fast precession so one can choose as the small time r the (approximate) 
precession period around the magnetic field, r = . The condition becomes 

(25) 




FIG. 3: Spin splitting of energy spectrum in the Rashba model. The green arrows represent the direction of the momentum 
and the red arrows represent the direction of the spin. 



Consider the case in which a small damping term exists a(^f x s). The equations of motion for the three spin 
components are 

cls^ r) ( dSy ds z . 

~~dt ~ ~ s z ll y + a \~u~ Sz m~ s v> 



ds v n . ( ds z ds 
- = s z \l x + a[——s x - 
y dt 

ds f~\ f~\ { ds y . 

— s X 4*y Sy\l x -(- oty Sy — s x ) . 



,, — s z Sl x + a(— -r~s x r^s z ) (26) 

dt y dt dt ' 



dt y dt 

It is clear that all the changes in s x are at least of second order in small quantities, so s x can be treated as a constant. 
The solution is approximately 

s y {t) = jw(Mxt - e~ t/T sina x t) 

(27) 

s z (t) = — f [1 - e - t/T {cosnt- — sin fit)]. 



fl 2 " v ft 



T 



where r = a Q s ■ After the oscillations in s y , s z are damped, what remains is 



s y (t » r) = — t 

J (28) 
s z (t»r) = ^. 

The spin, which was originally in the plane, acquires a small steady-state out-of-plane component. This component 
depends on s x , and s x has the opposite sign on different sides of the Fermi surface. As a result, s z up and down spins 
travel in opposite directions and a spin-Hall current is established. This derivation, illustrated by Fig. 3, is another 
way to clarify the role of spin precession in steady-state band structure spin currents. 



IV. EXTRINSIC SPIN DENSITIES AND CURRENTS 



Until now the scattering potential has been treated as a scalar. However, in general the spin-orbit interaction makes 
a contribution to the disorder potential. The disorder potential including the scalar and spin-orbit parts takes the 
form 

U kk , = [l + i\v(kxk')]^. (29) 
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FIG. 4: An electric field displaces the Fermi surface and the electron spins are tilted up for p y > and down for p y < 0. 



with lAkk' the matrix element of the scalar part of the potential. The second term comes from spin-orbit coupling 
in Ukk> and A is a material-specific constant. Configurational averages of terms of the form U^FU^' ll U^?l are also 
linear in rii but are rather lengthy to be displayed explicitly. 

The spin-orbit interaction also produces a correction 92] to the position operator r = f rd + 2A er x k where r or d 
is the ordinary position operator. Since the external electric field enters through the Hamiltonian He — eE ■ r the 
correction to r gives a side jump term H E = 2e\cr -kxE and this side jump term makes an additional contribution 
to the velocity operator v J E — —^jf- cr x E. In the kinetic equation extra spin- dependent driving term —(i/h) [H J E , fok] 
must be taken into account due to the change in the position operator. This term is nonzero if the equilibrium density 
matrix fok is spin-dependent. The side jump mechanism is typically classified as extrinsic although its contributions 
are manifold. These mechanisms were studied extensively in the anomalous Hall effect [111, H3, foiL [921 and recently 
their role in steady-state spin densities and currents have been highlighted in semiconductors [73, [ll, [79|, HU and 
metals. 01 

The solution for is found as above and the spin-dependent scattering termJ s acts on Ueh and produces an 
additional source terms for Seh, such that the RHS of Eq. (|20f becomes 

1 

S B fc - J s {n E k) - t [H E , S oh \ - Jj(foh)- (3°) 

Skew scattering emerges in the second Born approximation as a third-order term in the potential 14. Substituting for 
nEk and introducing the angles 71 = 9' — 8, 72 = 6" — 9' , 73 = 9 — 9" and the solid angles u/, oj" , 

f 6TT 2 \n t eT p m* 2 k 2d ~ 2 ~ 
Js\nk) = ^ o(k — Kf) cr ■ k x I 

(31) 



The integral over solid angles I is independent of k. In two dimensions the integrand is expanded in Fourier harmonics 
and the integration is straightforward. In three dimensions the integrand is expanded in Legendre polynomials of 71, 
72 and 73, and the integral over the two solid angles of Pi (71 )P m (72)^(73) is independent of k. The modification of 
the position operator produces an additional side jump scattering term 

-J j (fo k ) = - Hi6{£k ~ £F) . (32) 

Tp 

This term is related quite literally to transverse jumps undergone by a carrier during scattering [92| and its form 
reflects the conservation of momentum during such a side jump scattering event. [79|, [92j 
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The intrinsic source T, E k was discussed in the previous section. The extrinsic source is henceforth denoted by 
T Ek = — i [H 3 E , Sok] — Jj(fok) ~ Jsi^Ek)- If the band structure spin-orbit interactions are zero one finds easily 

Seh = h e S(e k - e F ) - J s {n Ek ) t p , (33) 

where the superscript ext indicates that only extrinsic contributions are considered. This expression averages to zero 
over directions in momentum space and does not give a spin density. It does however contribute to the spin current, 
giving a spin-Hall conductivity oc n^ 1 due to skew scattering and a spin-Hall conductivity ne\ independent of m due 
to side jump. These terms were discussed by Engel et al. [Z3] and by Hankiewicz and Vignale (79j |. 

If band structure spin-orbit interactions are nonzero then, as was pointed out by several groups, [8(| HH H2, US] 
spin precession is crucial in establishing the steady state. To illustrate this consider once again the decomposition 
^Ek = ^Ek\\ + &Ek± ™ to nnear ly independent components. In two-dimensional systems grown along (001) the entire 
extrinsic source term is orthogonal to H so , therefore Tehw — an d there is no term in that is oc n" 1 . The 
solution is 

next _ <T -n k X T Ek 

SEk± ~ 2^ • [S } 

The contributions to S^S^ due to skew scattering and side jump contain k to an even power, therefore skew scattering 
and side jump do not contribute to the spin current in the presence of spin precession (they can be restored by an 
external magnetic field B \\ z as found in Ref. [[80]].) The only contribution to the spin current comes from the 
spin-dependent driving term. The spin-Hall conductivity originating from this term is neX regardless of the form of 
the band structure spin-orbit interaction. Skew scattering in the presence of band structure spin-orbit interactions 
does give a steady-state spin density, which was pointed out in [ 83]], and a similar spin density arises from side jump. 

It is therefore seen that skew scattering and side jump contribute very differently to the steady-state density matrix 
when band structure spin-orbit interactions are present. This can be understood from the following argument. In 
order to get a sizable spin current from skew scattering, spins must be conserved as they travel. Scattering processes 
produce a separation between spin-up and spin-down, while precession tries to destroy the orientation of the spins. 
Thus once the spins are scattered it is the conserved spin fraction that carries the spin current. Evidently the same 
argument applies for side jumps undergone during scattering. 



V. DEFINITIONS OF THE SPIN CURRENT AND BOUNDARY CONDITIONS 

The spin current is defined in an intuitive manner by 3° — (1/2) {s 17 , v 1 }, and this definition is used by the majority 
of researchers working on this topic. Nevertheless, in the presence of band structure spin-orbit interactions this spin 
current is not conserved. The equation of continuity satisfied by the spin density and current was shown by Shi et al. 
to be different [32| from the usual equation of continuity for the charge density and current. A source term exists in 
this equation which reflects spin non-conservation. As a result, Shi et al. as well as Bryksin and Kleinert [34| have 
proposed an alternative definition of the spin current according to which 3° = d/ dt {fiS a ). The equation of continuity 
satisfied by this current still contains a source term, but this source term vanishes in many commonly used models. 
This definition was used in Refs. [[H, [69j]]. Recently, Tokatly [H[ has argued that the intuitive definition of the spin 
current represents a dissipative current which is conjugate to an effective SU(2) electric field. 

In order to obtain the spin accumulation the kinetic equation needs to be supplemented by boundary conditions. 
This was done by a number of groups. [H, [H, E3, EH E2, HE EH EH Unfortunately, it was shown that the form of 
the spin accumulation depends on the theoretical model of the boundary. [43], E3 The nontrivial physics associated 
with boundary conditions and the implications of various formulations of boundary conditions for spin transport were 
discussed by Bleibaum [39| and Tserkovnyak et al. [38[ An interesting argument due to Adagideli and Bauer [45[ points 
out that, in systems in which the band structure spin current is expected to vanish in the bulk (i.e. for Rashba and 
Dresselhaus spin-orbit coupling), the situation is not the same near the edges and near metal contacts, where a finite 
spin current exists which can be detected. 



VI. CRYSTAL SYMMETRY 



Spin densities and spin currents in a crystal are closely tied to the symmetry of the underlying lattice. An analysis 
relating the response tensor to the symmetry of the underlying crystal lattice has been enlightening in the context of 
nonequilibrium spin densities excited by an electric field. If the response of the spin density s to an electric field E is 
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given by s a = Q?Ej, nonzero components for the material-specific spin density response tensor are permitted only 
in gyrotropic crystals ||. For spin transport such an analysis was performed in Ref.[[36|]], determining the components 
of the spin-current response tensor allowed by symmetry in an electric field and providing systematic proof that spin 
currents in response to an electric field can be much more complex than the spin-Hall effect [4£j. This result is 
completely general and is not sensitive to the definition of the spin current or to whether the electric field is constant 
or time-dependent. 

The spin current operator can be defined as J? = ^(s a Vi + Vis' 7 ) or J? = d/dt (fiS a ). From a symmetry point of 

view these two definitions are equivalent. The spin current 3 is a second rank tensor that can be decomposed into a 
pseudoscalar part, an antisymmetric (spin-Hall) part, and a symmetric part. The pseudoscalar part is tr( jj) — | s ■ v 
and represents a spin flowing in the direction in which it is oriented. The symmetric and antisymmetric parts are 
given, respectively, by ^(s a Vi ± v a s l ). The pseudoscalar and symmetric parts will be referred to as non-spin-Hall 

currents. Under the full orthogonal group only the antisymmetric (spin-Hall) components of 3 are allowed, indicating 
that these components are always permitted by symmetry. 

In general, the spin current response of a crystal to an electric field E is charac teriz ed by a material tensor T defined 
by J° — TfjEj. For the 32 crystallograph y poi nt groups the symmetry analysis [l00j | for the tensor T is established by 
means of standard compatibility relations [10l| . One is particularly interested in those groups in which non-spin-Hall 
components may be present. The pseudoscalar part of the spin current is only allowed by 13 point groups, while the 
symmetric part is allowed in all systems except those with point groups O, Td (zinc blende), or Oh (diamond). Lower 
symmetries, allowing non-spin-Hall currents, are characteristic of systems of reduced dimensionality. In such systems 
it was shown that non-spin Hall currents exist. 36] For a quantum well grown along (113), the effective magnetic field 
corresponding to the Rashba spin-orbit interaction has a component in the ^-direction as well, and Ref. [[361]] showed 
that the spin conductivities a xx , a yy , a yx and a xy are nonzero. 

Skew scattering and side jump contributions to the spin current in the presence of band structure spin-orbit 
interactions can also be restored by growing the structure along a direction that is not one of the main crystal axes. 
The contribution to the spin density matrix due to extrinsic mechanisms has two parts, 

next 1 ( T~EkzQkz\ ,„_,. 

Bfc|1 2 V — & Jv^kT (35) 

and S^kX' which ls given by Eq. (|3"4"|) . Taking again a quantum well grown along (113), Sgj£n gives no spin density 
but gives a spin current from skew scattering and side jump. Skew scattering and side jump give the same nonzero 
components of the spin conductivity tensor, a xx , a yy , a yx and a xy . 



VII. EXPERIMENTAL SITUATION 



Spin densities in semiconductors were first predicted by Ivchenko [4[ and observed the following year in tellurium 
by Vorob'ev. || Recent experiments [l^, [HI El have also reported the observation of a steady-state spin density in 
semiconductors. Experimentally the spin density can be found by Kerr rotation, in which a beam of linearly-polarized 
light is sent into the sample and its polarization vector rotates by an angle that is proportional to the spin polarization. 

In so far as measuring a spin current, the situation is more complicated. Initially the approach adopted was to 
generate a spin current which would flow to the edges of the sample where it would produce a spin accumulation. 
This spin density at the edge of the sample could then be detected by Kerr rotation or an equivalent technique. 
J. Wunderlich et al. [H[ used photoluminescence to detect a spin accumulation due to a spin-Hall current in a two 
dimensional hole gas, which is believed to be due to spin-orbit interactions in the band structure. Kato et al. (l9j 
used Kerr rotation to detect a spin accumulation due to a spin-Hall current in rt-GaAs, in which the spin current is 
believed to be due to extrinsic mechanisms [77| . A similar experiment was carried out shortly thereafter by Sih et al. 
[20| . Stern et al. [2l[ used the same technique to detect a spin-Hall current in ZnSe at room temperature. Following 
this work, Sih et al. [22| performed an experiment designed to demonstrate explicitly that the observed edge spin 
accumulation was due to the spin-Hall effect, as opposed to edge effects associated with the electric field. The group 
manufactured samples with a series of transverse channels, such that the spin-Hall current generated by the electric 
field was allowed to drift into regions where the electric field was effectively zero. Since a spin accumulation was still 
measured at the edge, this showed unambiguously that the effect was due to the spin current. Further experiments 
by Stern et al. [22j| imaged the spatial distribution of the spin accumulation generated by the spin-Hall effect as well 
as its behavior in a magnetic field. Chang et aZ.(2~i| also reported, using photoluminescence, a spin accumulation due 
to the spin-Hall effect in InGaN/GaN superlattices. 

A second type of experiments relies on an effect referred to as the inverse spin-Hall effect. Briefly, a spin current in 
turn generates a charge current, which can be detected by conventional means. This was shown by Hirsch [l4| and its 
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extension to a Landauer-Buttiker type multi-terminal nonlocal measurement was presented by Hankiewicz et al. [74j 
This technique was used by Valenzuela and Tinkham p5| to detect a charge current as a result of the spin current 
in aluminium, and by Saitoh et al. [13] in platinum. Kimura et al. [28[ observed the spin-Hall effect in platinum at 
room temperature and their findings were explained theoretically by Guo et al., (76j who demonstrated that the effect 
is due to band structure spin-orbit interactions. Vila et al. also obtained results for platinum nanowires [29j while 
Scki et al. (30l | reported a giant spin Hall effect in FePt/Au devices. The group used a multi-terminal device with a 
gold cross and FePt acting as a spin injector. The spin-Hall resistance was measured to be 2.9 mfi and is attributed 
to the large skew scattering in gold, being thus extrinsic in nature. Weng et al. [3l| carried out similar experiments 
on platinum, aluminium and gold and obtained results in agreement with those found to date. It should be noted 
that Cui et aL(26| also observed an electrical current induced by an optically-generated spin current. 



VIII. FUTURE DIRECTIONS 



Whereas the scientific community working on steady-state spin densities and currents appears to be in agreement 
that spin currents exist and are experimentally measurable, a number of questions remain to be addressed in the 
future. For example, the relative magnitude of band structure spin currents and spin currents due to extrinsic 
mechanisms such as skew scattering and side jump remains to be determined for a general band structure spin- 
orbit interaction. In addition, the fact that no unique definition of the spin current exists causes difficulties in the 
comparison of experimental data with theoretical predictions. This ambiguity is exacerbated by the fact that different 



definitions of the spin current give results that often differ by a sign. [32j, |33| Thanks to the non-conservation of 
spin, the relationship between spin current and spin accumulation at the boundary remains to be clarified. It appears 
that what happens at the boundary is sensitive to the type of boundary conditions assumed. [H, Q Thus so far as 
quantitative interpretation of experimental data is concerned, theory still has some way to go. 

On the experimental side, despite tremendous progress, the community is still searching for a reliable way to 
measure, as opposed to detect, spin currents directly. A possible new experimental avenue relies on magnetoresistance 
caused by an edge spin accumulation, proposed by Dyakonov. [85] A spin current produces a spin accumulation near 
the sample edges, which in turn causes the sample resistance to decrease by a small amount, whereas an external 
magnetic field can destroy the edge spin polarization and yield a positive magnetoresistance. An alternative path was 
followed by Wang et aL,[46( who started from the Dirac equation and obtained in the weakly relativistic limit a set of 
Maxwell equations in the presence of spin-orbit interactions. Although the spin current does not appear explicitly, it 
is contained in the Maxwell equations and the authors demonstrate that the relativistic conservation laws imply that 
the spin current yields an electrical polarization, which could be detected directly. Similarly, a set of Maxwell-like 
equations for spin was formulated by Bernevig et a£.[47j These methods remains to be attempted experimentally. 

Practically, the question of what to do with the electron spin once it has been generated or transported to the 
edge of the sample remains. The revolutionary electronic device that harnesses spin currents for a practical purpose 
remains to be made, and the challenge of its design confronts experimentalists and theorists alike. 

The research at Argonne National Laboratory was supported by the US Department of Energy, Office of Science, 
Office of Basic Energy Sciences, under Contract No. DE-AC02-06CH11357. 
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